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ABSTRACT 



> 
CO 

in 

■ Context. B[e] supergiants are known to have non-spherical winds, and the existence of disks that are neutral in hydrogen close to their stellar 
f^*D surface has been postulated. A suitable mechanism to produce non-spherical winds seems to be rapid rotation, and at least for three B[e] 
\Q . supergiants in the Magellanic Clouds rotation velocities at a substantial fraction of their critical velocity have been found. 

' Aims. We want to find suitable recombination distances in the equatorial plane of rapidly rotating stars that explain the observed huge amounts 
"p^ \ of neutral material in the vicinity of B[e] supergiants. 
Q . Methods. We perform ionization structure calculations in the equatorial plane around rapidly rotating luminous supergiants. The restriction 
to the equatorial plane allows us to treat the ionization balance equations 1-dimensionally, while the stellar radiation field is calculated 
2-dimensionally, taking into account the latitudinal variation of the stellar surface parameters. The stellar parameters used correspond to those 
known for B[e] supergiants. The assumptions made in the computations all have in common that the total number of available ionizing photons 
at any location within the equatorial plane is overestimated, resulting in upper limits for the recombination distances. 

Results. We find that despite the drop in equatorial surface density of rapidly rotating stars (neglecting effects like bi-stability and/or wind 
compression), hydrogen and helium recombine at or close to the stellar surface, for mass loss rates M 2; 5 x lCr 5 M G yr~' and rotation speeds in 
excess of Vrot,eq/vcrit - 0.8. 



o 

CO 



X 



Key words. Stars: rotation - Stars: mass-loss - Stars: winds, outflows - supergiants 



1. Introduction 

In a series of papers Maeder & Meynet discussed the impor- 
tance and the influence of rapid rotation on the evolution of 
stars, the chemical yields and non-spherical mass and angu- 
lar momentum loss (see e.g Maeder, 1999, Meynet & Maeder, 
2000 Maeder & Meynet, 2000|>. The influence of rotation on so 
many stellar parameters results also in the shaping of the wind 
and the nebula. There has been the suggestion that the appear- 
ance of non-spherical winds around some massive and lumi- 
nous stars might be caused by rotation (Maeder 2002; Maeder 
& Desjacques 2001 1. Stars for which this might be appropriate 
are the Luminous Blue Variables and the B[e] supergiants. 

Zickgraf et al. £T985l suggested that the hybrid charac- 
ter of the optical spectra of B[e] supergiants is due to a non- 
spherical (two-component) wind. The strong infrared excess 
indicates the presence of a huge amount of hot circumstellar 
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dust, and polarimetric observations, e.g. by Magalhaes ( 1992 1, 
Magalhaes et al. (2006 1, and Melgarejo et al. J2001> . con- 
firmed the non-spherical geometry of the circumstellar mate- 
rial around B[e] supergiants. The existence of a geometrically 
thick circumstellar disk responsible for the polarized emission 
and the location of the hot dust seems nowadays to be well es- 
tablished (for an overview see e.g. Kraus & Miroshnichenko 
2006 1. The formation mechanism of these disks is, however, 
still rather unclear. There exist two promising approaches: (1) 
the bi-stability mechanism introduced by Lamers & Pauldrach 
( 1991 ) and further investigated, especially with respect to the 
influence of rotation on the formation of B[e] supergiant stars' 
disks, by Pelupessy et al. (2000 1, and (2) the wind-compressed 
disk, introduced by Bjorkman & Cassinelli dl993> . Both mod- 
els however have difficulties in explaining all observed quanti- 
ties of the B[e] supergiants' disks in a self-consistent way. In 
addition, there is still no consensus about the nature of these 
disks, whether they can be described by a high density equa- 
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torial outflowing wind or by a Keplerian viscous disk (see e.g. 
Porter l2UU51 Kraus & Miroshnichenko l2006V 

Recently, these disks have been suggested to be neutral in 
hydrogen in the vicinity of the stellar surface. Tracers for hy- 
drogen neutral material are e.g. the strong [Oi] emission lines 
arising in the optical spectra of B[e] supergiants. Modeling of 
their line luminosities revealed that, in order to keep the mass 
loss rate of the star at a reliable value, these lines must origi- 
nate within a few stellar radii from the surface (Kraus & Borges 
Fernandes 2005, Kraus et al. 2006 1. In addition, several B[e] 
supergiants are found to show band-head emission from hot 
(3000 - 5000 K) CO gas (McGregor et al. IT9881 . Follow-up 
studies of high-resolution spectra for at least one of them led to 
the conclusion that this hot CO gas is located at about 2-3 AU 
from the hot stellar surface (Kraus 1277001 Kraus et al. 2000 1. 
The existence of neutral material close to these luminous B[e] 
supergiants is surprising and needs to be investigated in detail. 
The goal of our study is therefore to find scenarios that allow 
neutral material to exist close to the surface of these stars. 

In a first attempt, Kraus & Lamers ( 2003 , hereafter Paper 1) 
calculated the ionization structure of B[e] supergiants, assum- 
ing a latitude-dependent mass flux that increases from pole to 
equator. With such a model they could show that, even with 
moderate total mass loss rates, hydrogen recombines in the 
equatorial direction close to the star leading to a hydrogen neu- 
tral disk-like structure. Here, we investigate the influence of 
rotation on the stellar parameters and consequently on the ion- 
ization structure in the winds and disks of B[e] supergiants. 

Rotation causes a flattening of the stellar surface and there- 
fore a reduction of the local net gravity in the equatorial re- 
gion. The decrease in gravity from pole to equator equally re- 
sults in a decrease of the stellar flux which is proportional to 
the local net gravity. Hence, the effective temperature also de- 
creases from pole to equator, known as gravity darkening (or 
polar brightening, von Zeipel lT924> . The latitude dependence 
of the gravity and effective temperature has also impacts on the 
stellar wind parameters (see e.g. Lamers & Cassinelli ll999> : 
The escape velocity following from the balance between gravi- 
tational and centripetal forces becomes latitude-dependent, de- 
creasing from pole to equator. The same holds for the terminal 
wind velocity which is (for line-driven winds) proportional to 
the escape velocity. Even the mass flux from the star tends to 
decrease from pole to equator if gravity darkening is taken into 
account in the CAK theory as shown by Owocki et al. (09981. 
More important for the ionization structure calculations is the 
density in the wind, and we will show in Sect.[2]that the density 
at any given distance also decreases from pole to equator. A 
rotating star will therefore have a less dense wind in the equa- 
torial region, unless special effects such as bi-stability or wind 
compression play a role. 

Both important parameters in the ionization balance equa- 
tions, i.e. the effective surface temperature and the surface den- 
sity, decrease from pole to equator. While the decrease in sur- 
face temperature tends to decrease the number of available pho- 
tons suitable to ionize H and He, the decrease in surface density 
reduces the optical depth along the line of sight from a point in 
the wind to the star. Both effects are therefore counteracting 
with respect to the location where recombination takes place. 



While a reduction of ionizing photons will shift the recombi- 
nation distance towards the star, the reduction in optical depth 
along the direction to the star will shift it further outwards. The 
outcome of the ionization balance equations is therefore unpre- 
dictible and very sensitive to the chosen input parameters. We 
thus investigate the ionization structure in the wind of a rotating 
star in more detail. 

The paper is structured as follows: In Sect.|2]we provide 
the equations that describe the surface distribution of the effec- 
tive temperature, mass flux, escape velocity (and hence termi- 
nal wind velocity), and hydrogen density for a rigidly rotating 
star. The ionization structure calculations restricted to the equa- 
torial plane of the systems are performed in Sect.|3]where also 
the results for the recombination distances of helium and hy- 
drogen are shown. The influence of the assumptions on these 
results and the applicability of the models to B[e] supergiants 
are discussed in Sects. |4] and |5] respectively, and the conclu- 
sions are given in Sect.|6] 

2. The surface and wind structure of rigidly 
rotating stars 

In this paper, we restrict our investigations to rigid rotation 
only, and we neglect any influences due to bi-stability and 
wind-compression. 

2.1. The shape of the stellar surface 

The potential <t> of a rotating star is given by the sum of gravita- 
tional and centrifugal potential. The latitude dependence of the 
stellar radius, R{9), i.e. the shape of the star, is determined by 
the equipotential surfaces, <$>{R{9), 9, (/>), for which it is assumed 
that all the mass is concentrated in the core. These equipotential 
surfaces are given by 



®(R(6), 6,<f>) = - 



GM eff _ 1 

R(9) 2 



R 1 (9)Q. 1 sin 2 (6>) 



(1) 



where 9 is the co-latitude with 9 — at the pole, M s g is the ef- 
fective stellar mass, i.e. the stellar mass reduced by the effects 
of radiation pressure due to electron scattering, and £2 is the an- 



gular velocity. With the definitions of x(9) 



R(9)/R< 



eq> 



-^(GM e ff)//?eq, and to = v ro t,eq/v c rit, the latitude-dependent stel- 
lar radius is found from Eq. {0 which results in the following 
cubic function 
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x(9 + 0) 
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(3) 



(4) 



Equations Q and @ describe the stellar radius at all latitudes 
for a star rotating rigidly with a specific value of u>. 
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2.2. The latitude-dependent surface temperature 

Rotation not only influences the radius of the star but also re- 
sults in a latitude-dependent surface temperature distribution, 
because the stellar flux, F, is proportional to the local effective 
gravity, g e g, which is calculated from 



?eff 



= -V«D 



1 



R 2 (0) 



1 - 



— j^u? sin 2 

^eq 



(5) 



Since F(G) = o-T* s (6) 
behaves as 



g e ff , the surface temperature T e g (6) 



Tim 



i 



R 2 (ff) 



1 - 



2 ■ 2 , 

a) sin ( 



(6) 



or, if we express the latitude-dependent effective temperature 
in terms of the polar temperature, r e ff(pole), replace R(9) by 
x(0), and make use of relation 0, 



r e 4 ff (pole) 



= ^(pole) 



fl 2 (pole) 

~rW 



i 



R (0) 2 . 2 

— sin 2 9 
Rl, 



1 (1 - x 3 (6)uj 2 sin 2 6) 
xHff) (1 + I w 2 ) 2 



(7) 



For the purpose of our paper it is important to treat the sur- 
face effective temperature (as well as all other following pa- 
rameters) properly, which means that we have to take into ac- 
count the rotationally distorted stellar surface. How different 
the results can be when accounting only for gravity darken- 
ing but neglecting the real shape of the star is shown in Fig.^ 
There we compare the surface temperature distribution calcu- 
lated from Eq. which accounts for the distorted surface with 
the one resulting from a rotating star but under the assumption 
of an unperturbed, spherical surface (i.e. R(0) = R eq = R). In 
this latter case, T e ff(9) resulting from Eq. (|6} is simply given by 
(see e.g. Lamers & Cassinelli 1999, Lamers 1 120041 



T 4 eS (0) = T 4 eff (pol S )(l-cj 2 sm 2 o) . 



(8) 



This equation also describes globally the influence of rigid 
rotation, i.e. the drop in temperature from pole to equator. 
However, the absolute value of the effective temperature at any 
location on the stellar surface is different, as is obvious from 
the comparison in Fig. [2 

- For to £ 0.8 the rotationally distorted surface is cooler at 
all latitudes. 

- For oj ^ 0.8 the effective temperature of the rotationally 
distorted surfaces is higher for small to intermediate lati- 
tudes, but becomes (much) lower in the equatorial regions, 
compared to the corresponding spherical surfaces. 

These severe differences in surface temperature distribution 
have non-negligible effects on the stellar radiation field at any 
point in the wind. A proper treatment of the stellar parameters 
by accounting for the rotationally distorted surface is therefore 
an important ingredient in our ionization balance calculations. 



1 Please note the typo in Lamers' Eq. (5) where it should be cos 2 (0) 
instead of cos(£>), and that in his paper 6 is measured from the equator. 
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Fig. 1. Effective temperature distribution on the surface of a 
rotating star. The different curves, which are normalized to 
the polar temperature, are for different rotational velocities in- 
dicated by co from low values (upper curves) to high values 
(lower curves). In the top panel the rotational distortion of the 
stellar surface has been neglected. 

2.3. The latitude-dependent mass flux 

For the mass flux, F m we follow the description of Owocki et 
al. ( 119951 given by their Eq. (2) 



F m (6) 
F m (pole) 



F(ff) 



F(pole) 



geff(pole) 



(9) 



This equation describes the latitude-dependent mass flux ac- 
cording to CAK theory (Castor et al. H975i . Neglecting bi- 
stability effects, which means that the force multiplier a is con- 
stant all over the surface, and introducing gravitational darken- 
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ing according to the von Zeipel theorem (i.e. F(0) 
results in 



<eff 



(0) 



1 (1 - x 3 (6)u) 2 sin 2 6) 



F m (ff) = = 

F ra (pole) geff(pole) x 2 (0) (1 + ±« 2 ) 2 



(10) 



The surface distribution of the mass flux for different values of 
a) is shown in the top panel of Fig.|5] 

2.4. The latitude-dependent terminal velocity 

The escape velocity of a rotating star follows from balancing 
gravitational and centripetal forces on the stellar surface which 
means that the effective gravity g e s(0) must equal v 2 sc (#)//?(0). 
Since the terminal wind velocity, v^, is, according to line- 
driven wind theories (see e.g. Lamers & Cassinelli 119991 
Chapter 8), proportional to the escape velocity, v esc , we find the 
following relation for the latitude dependence of the terminal 
velocity 



Voo(f) = v ro (pole) 



(l - x\0)cj 2 sin 2 6) 
(x(0)(l + i« 2 )) 1/2 



1/2 



(ID 



The latitude dependence of the terminal velocity for different 
values of o> is shown in the middle panel of Fig. [2 

2.5. The latitude-dependent wind density distribution 

It is known that the mass flux and the terminal velocity de- 
crease from pole to equator, even if the rotational distortion of 
the stellar surface is neglected in their derivation (see Lamers 
& Cassinelli 1999 1. For the ionization structure calculations, 
however, we need to know the density distribution in the wind. 

In a non-rotating, spherically symmetric stationary wind, 
the density at any location r in the wind is related to the mass 
loss rate, M, of the star and the wind velocity, v(r), via the 
equation of mass continuity 



nn(r) 



M 



R 1 



47Tfminr 2 v(r) fimnv(r) r 2 



(12) 



where /j is the mean molecular weight and «h denotes the parti- 
cle density of hydrogen given in cirT 3 . From the right-hand side 
of this equation it follows immediately that in a non-spherically 
symmetric wind, the radial density distribution at any latitude 
can be written in the form 



F m (0) R 2 (9) 



fimnv(6, r) r 2 



(13) 



For our further calculations, we assume that the wind veloc- 
ity is constant in radial direction, i.e. v(0,r) = v{0,R) = v(ff). 
Therefore we can re-write Eq. ill 3b in the form 



n H (0,r) = n H (6,R(ff)) 



R 2 (0) 
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Fig. 2. Distribution of the mass flux (top panel), the escape ve- 
locity (mid panel), and the hydrogen density (bottom panel) on 
the surface of a rotating star. The different curves in each plot, 
which are normalized to the corresponding polar value of the 
parameters, are for different rotational velocities indicated by 
u>. All parameters drop from pole to equator. This effect be- 
comes stronger with increasing values of u. 



where 



nn(0,R(ff» 



F m (0) 
um H v(6) 



(15) 



defines the density distribution along the stellar surface. 

With the additional simplification of v(0) = v 00 (6 l ), we can 
express the surface density distribution, rtu(0, R(9) ), by using 
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the relations for the mass flux given by Eq. and for the 
terminal velocity given by Eq. (1111 



nn{e,R{6)) 



F m (pole) (l -x 3 (0)w 2 sin 2 0) 



1/2 



^m H v M (pole) + 

(l -x 3 (0)co 2 sin 2 6») 1/2 



3/2 



n H (pole) ■ 



(x(0)(l + \u> 2 )) 



Hi 



(16) 



(17) 



This surface density distribution, as the result of the ratio of 
mass flux to terminal velocity, is plotted for different values 
of to in the lower panel of Fig.|2] It also decreases from pole 
to equator. This means that a rigidly rotating star will have a 
less dense wind in the equatorial region, unless bi-stability and 
wind compression play a role. 

3. Ionization structure calculations 

Since we are searching for the existence of a hydrogen neutral 
equatorial region, we restrict our calculations to the equato- 
rial plane only. This leads to the simplification of a symmetri- 
cal stellar radiation field with respect to the equatorial plane. 
Therefore, it is sufficient to solve the ionization balance equa- 
tions along one radial direction, which we will call the y-axis. 

As in Paper 1 , our model wind consists of hydrogen and 
helium, only. This means that we have to solve two coupled 
ionization balance equations that are treated in the on-the- 
spot (OTS) approximation. This approximation states that ev- 
ery photon generated via recombination and able to ionize hy- 
drogen or helium will be absorbed immediately in the close 
vicinity of its generation location. The ionization balance equa- 
tions for this case are given in Sect. 4. of Paper 1 . The recom- 
bination distance is found by applying a root-finding routine. 
Usually, a few iteration steps are sufficient for an accuracy in 
distance better than 1%. 

3.1. The stellar radiation field 

Besides the OTS approximation, which defines the diffuse ra- 
diation field, we need to calculate the stellar radiation field at 
any point in the equatorial plane, or, due to the symmetry in 
our case, at any point along the y-axis. Differently from the 
treatment in Paper 1 we no longer use the assumption that the 
star is a point source. The stellar parameters of a rotating star 
vary strongly over the stellar surface, especially for increasing 
stellar rotation. We therefore calculate the stellar radiation by 
integrating the latitude-dependent surface flux over the rota- 
tionally distorted stellar surface. This is done in the following 
way: 

- We define the stellar input parameters T e ff (pole), g e ff (pole), 
Voo(pole), /?*(sphere), F m (pole). 

- We define the rotation velocity, to. 

- With co and /^(sphere) we calculate the shape of the stellar 
surface, i.e. x(0). We thereby make use of Eq. (0 and of the 
mass conservation that relates the spherical radius to the 
equatorial and polar radii via = /? 2 q ^ po i e . 



- At each location r along the y-axis we determine the angu- 
lar extent of the stellar surface and its shape. This defines 
the size of the stellar surface (or the surface segment) from 
which radiation will arrive at point r. 

- Along this stellar surface segment we calculate the distri- 
bution of T e ff(ff) and g e ff (&), and the resulting radiation tem- 
perature T m d{&) (see Sect. l3.2> . 

- The stellar flux as a function of latitude is then approxi- 
mated by B v (T rdd (e)). 

- The total stellar radiation field at point r along the y-axis 
follows from integration of B v (T rd d(0)) over the segment of 
the rotationally distorted surface. 

- To account for optical depth effects, we calculate the mini- 
mum optical depth, which occurs along the y-axis (because 
of the shortest distance and the lowest density). This op- 
tical depth is adopted for all directions towards the stellar 
surface. 

In this calculation of the stellar radiation field at any point r 
along the y-axis we make one important approximation, which 
is the adoption of the minimum optical depth towards all di- 
rections from r to the stellar surface. This assumption results 
in an overprediction of the available ionizing photons, because 
the stellar radiation from the higher latitudes with higher J e ff 
will be absorbed less. The resulting recombination distance 
will therefore be overestimated. 



3.2. The radiation temperature 

The stellar radiation temperature is defined in our calculations 
as the Planck temperature which describes the part of the stellar 
spectrum that delivers the ionizing photons, i.e. the spectrum 
shortwards of 912 A. The determination of the latitude depen- 
dence of the radiation temperature is not straightforward. We 
therefore briefly explain how we calculate it. 

We start with Kurucz model atmospheres (Kurucz fl979i 
for solar metallicity stars and log g values between 2.0 and 3.5. 
We do not investigate higher values of log g because we are 
mainly interested in giants and supergiants. For each model 
atmosphere in this logg range and for all available effective 
temperatures we fitted a Planck function to the spectrum short- 
wards of 912 A. This delivers the radiation temperature for the 
corresponding (r e ff, g e g) combination. The radiation tempera- 
ture is therefore a function of these two parameters, i.e. r ra( j = 
7rad(7eff,geffX an d we found the following useful parametriza- 
tion: 

log T mi = A(logg) log r eff + BQogg) (18) 

with the two functions A(logg) and fi(logg) given by 

A(logg) = 1.222 - 0.058 logQogg- 1.9) (19) 
flQogg) = 0.235 log(logg - 1.89) - 1.13 (20) 

in the range of 8 000 £ r eff < 30 000 K. The error in radiation 
temperature introduced by this fitting procedure is less than 5% 
for the higher values of T e ff, and less than 3% for the lower 
ones. We thus can compute T m d(6) for any combination of T e g 
and geff values that will occur along the surface of a rigidly 
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rotating star, and hence we can calculate the appropriate stellar 
radiation at all locations on the stellar surface. 

Since T m ^(ff) is a function of both T e ff and g e g , its latitude 
dependence will be different from the one of the effective tem- 
perature (see Sect. l3.4> . and will also vary for stars with differ- 
ent stellar parameters. 

3.3. Description of the model stars and their winds 

For our calculations we chose stars with a polar effective 
temperature of r effjPO i e ^ 24 500 K and log g eff , po i e = 3.5. 
According to Eq. (II 81 this combination results in a polar ra- 
diation temperature of T la ^ po \ e ^ 17 000 K. 

The radius of the non-rotating star is fixed at R, = 82 R . 
Together with the chosen effective temperature this results in a 
luminosity of the non-rotating star of L* =s 2.2 x 10 6 L , which 
places the star in the B-type supergiant region within the HR 
diagram. 

The polar values of the effective temperature and gravity 
(and hence radiation temperature) are fixed for all our model 
calculations. Fixing the polar effective temperature means that 
we are not calculating a star that is spinning up. This would re- 
sult in an increase of polar temperature with increasing rotation 
velocity. Instead, we are calculating stars with the same polar 
effective temperature having different rotation velocities. This 
means that we are dealing with stars of different luminosities. 
The total luminosity of a star is 



n 2n 



J J o-Tt 9 {6)dS 



(21) 



o o 



where dS is the surface element of an ellipsoid, given by 



dS 



cos 2 8 + R 2 , sin 

pole 



0) R eq sin 0ddd<f> . 



(22) 



The difference in luminosity of rotating stars having the same 
polar effective temperature is shown in the upper panel of Fig. [3] 
where we plotted the distribution of the stellar luminosities as a 
function of u. With increasing rotation speed, the stellar lumi- 
nosity drops. The difference is largest between the non-rotating 
and the critically rotating star, and is about a factor 2. 

Similarly, we can calculate the mass loss rates of our model 
stars. The mass loss rate follows from 



7T 2k 



F m (ff)dS 



(23) 



o o 



and the ratio of the mass loss rate over the mass loss rate of the 
spherical star is plotted as a function of oj in the lower panel of 
Fig. |3] The mass loss rate shows a difference of about a factor 
2 between the non-rotating and the critically rotating star. The 
behaviour of both the luminosity and the mass loss rate with 
lo is identical because both T* s {8) and F m {8) are proportional 
to g e ff{6) (see Eqs. and dlOt respectively). We can therefore 
conclude that an almost critically rotating star will have only 
half the total luminosity and half the total mass loss rate of its 
non-rotating counterpart, having both the same polar mass flux 
and effective temperature. 
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Fig. 3. Top panel: Stellar luminosity as a function of to for 
our model stars with r e ff jP oi e = 24 500 K and R* = 82 7? . 
Bottom panel: Mass loss rate as a function of o>, normalized 
to the spherical mass loss rate. Both parameters show a de- 
crease with increasing rotation speed. The difference between 
the non-rotating and the critically rotating star is about a factor 
2. 

In our calculations, the polar mass flux, f m po i e , is a free pa- 
rameter. Its value is varied over a large range to investigate the 
ionization structure of rotating stars and to find the recombina- 
tion distances (see Sect. 13.4^ 

We further use a distant independent wind velocity, which 
is set to the terminal velocity, i.e. v{6, r) = v oo (0). For the po- 
lar wind velocity we adopt Vo^poie = 2000 km s -1 . The radial 
electron temperature distribution in the wind is assumed to be 
constant, and we set T e (6, r) = T e (ff) = 0.8 T e ff(ff). The influ- 
ence of these assumptions and simplifications on the results are 
discussed in Sect.|4] 

3.4. Recombination in the equatorial plane 

We calculated the equatorial recombination distance of H and 
He for stars with a large range in polar mass fluxes, F m;PO i e . For 
each mass flux we considered rotation velocities u> covering 
the complete range from to 1. In Fig.@]we show the results of 
three representative models, calculated for F m po i e = 1.0xl0~ 5 , 
1.5 x 10~ 5 , and 2.0 x 10~ 5 g s"'cm" 2 . The recombination radii 
for helium (left panel) and hydrogen (right panel) are given in 
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Fig. 4. Distance r in the equatorial plane in terms of the equa- 
torial radius R sq (u)) at which recombination of Hen (left panel) 
and Hn (right panel) takes place for stars rotating with different 
velocities, indicated by a>. The curves are (from top to bottom) 
for polar mass fluxes of: F m (pole) = 1 .0 x 10~ 5 , 1 .5 x 10~ 5 and 
2.0 x lO^gs-'cm" 2 . 



units of the corresponding (^-dependent stellar equatorial ra- 
dius, R e q(a>). Helium recombines for all models already close 
to the stellar surface. The recombination distance is found to 
decrease steadily with increasing stellar rotation, reaching the 
stellar surface for co ^ 0.7. This means that for rapidly ro- 
tating stars helium is neutral at the stellar surface. For hydro- 
gen, the situation is different. For polar mass fluxes F mjPO i e S= 
1.5 x 10~ 5 g s^'cm -2 the recombination distances decrease with 
increasing a>, reach a minimum in the range a> - 0.70 . . .0.75, 
and increase again for higher rotation speeds. If the input polar 
mass flux is higher than 1 .5 x 10 _5 g s^'cm -2 , hydrogen shows 
the same trend as helium, i.e. the recombination distances de- 
crease steadily for increasing stellar rotation. They reach also 
the surface of the star, but for rotation speeds a> ^ 0.85. The 
models with F miPO i e = 1.5 x 10~ 5 g s _1 cm _2 seem to be the case 
in "transition". They show a kink at a> = 0.75 and a subsequent 
steep drop in recombination distance. 

What causes the minimum and especially the strong in- 
crease in hydrogen recombination distance for the lower mass 
flux models? This rather unexpected behaviour can be under- 
stood upon inspection of Fig.|5] There we plotted the variation 
of the individual equatorial surface parameters with a>, and of 
special interest are the radiation temperature and the particle 
density. Fig.|5]shows that with increasing rotation velocity, the 
density drops much quicker from pole to equator than the tem- 
perature. For recombination to take place right above the stel- 
lar surface, the number of ionizing photons has to be reduced. 
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Fig. 5. Surface parameters in the equatorial plane, normalized 
to their polar values, as functions of the rotational velocity. «h 
is the density at any distance. 

This can be done either by decreasing the radiation tempera- 
ture, or by increasing the equatorial surface density and hence 
the optical depth. Since the decrease in radiation temperature 
is determined by the rotation velocity (with a fixed input po- 
lar value) we can only increase the input polar mass flux to 
achieve a higher surface density for a given rotation velocity. 
A higher density also has the effect of triggering recombina- 
tion, this is however only a secondary effect while the blocking 
of the raidation field is the more important one. The density 
in the top model shown in Fig .0] is no longer high enough for 
stars with cj > 0.7 to absorb the ionizing photons provided by 
the still rather high radiation temperature. Therefore, ioniza- 
tion takes over again and shifts the recombination distance for 
higher rotation velocities further out. Whether recombination 
takes place close to the star therefore sensitively depends on 
the chosen input parameters of the rotating star. 

4. Discussion 

For our calculations we made a few assumptions, and we briely 
discuss their influence on the model results: 

The electron temperature. The winds of hot stars are 
known to start with an electron temperature of about 0.8 T e g 
at the stellar surface. Further out, they cool quickly (within 
a few stellar radii) and converge towards a more or less con- 
stant (terminal temperature) value (see e.g. Drew 1989 1. The 
adoption of a constant (in radial direction) and maximum (i.e. 
T s (0, r) = 0.8 T e ff(ff)) electron temperature reduces (or sup- 
presses) the total number of recombinations taking place, be- 
cause the recombination coefficient is small for high tempera- 
tures but increases with decreasing temperature (see Fig. 2 in 
Paper 1). Our assumption therefore inhibits recombination and 
shifts the equilibrium of the ionization balance in favour of the 
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ionization of the wind material. This means that we have over- 
estimated the recombination distance. 

The wind velocity. The velocity distribution in winds of 
hot stars can usually be approximated by a /3-law (see e.g. 
Lamers & Cassinelli 1999 1 which describes the increase in ve- 
locity from the small surface value to the terminal velocity. 
First, using the maximum (i.e. terminal) velocity instead of the 
more realistic yS-law increasing velocity distribution results in 
an underestimation of the density, especially at distances close 
to the star. Second, the chosen value of 2000 km s for the po- 
lar terminal velocity is rather high for a B-type supergiant and 
results equally in an underestimation of the density. Both as- 
sumptions therefore lead to an underestimation of the optical 
depth seen by the stellar radiation and an overestimation of the 
recombination distance. 

The optical depth. We calculate the optical depth properly 
only along the y-axis. Since the surface density of a rotating star 
drops from pole to equator and since the distance from the stel- 
lar surface to any point along the y-axis is the shortest one over 
which stellar photons can be absorbed, this optical depth value 
is the smallest. Adopting this minimum optical depth for all 
directions towards the stellar surface therefore underestimates 
the real optical depth and allows more ionizing photons to pen- 
etrate the wind material to larger distances. This leads to an 
overestimation of the recombination distance. 

The OTS approximation. In the OTS approximation it 
is assumed that every photon created via recombination and 
able to (re-)ionize hydrogen or helium will be absorbed in the 
close vicinity of its generation location, and none will escape 
from the wind. These leads to additional ionizing photons (bet- 
ter known as the diffuse radiation field) at any location in the 
wind. The OTS approximation, even if not fully applicable for 
the lower density regions, tends to overestimate the number of 
available photons everywhere in the wind and therefore favours 
the ionization of the wind material. 

All these assumptions and simplifications made in our com- 
putations and listed here show the tendency to overestimate 
the number of ionizing photons available at a certain location 
in the wind. Consequently, the recombination radius is shifted 
away from the star which means that our calculated distances 
at which the material is found to recombine are upper limits. 

5. Applicability to the B[e] supergiants 

In Sect. l3.4l we showed that the equatorial winds of rapidly ro- 
tating stars might be neutral in hydrogen right from the stel- 
lar surface, even though the density in the equatorial wind is 
much lower than in the polar regions. Since our main goal is to 
find possible formation mechanisms for hydrogen neutral disks 
around B[e] stars and especially B[e] supergiants, we discuss 
here how reliable the results are and whether they are indeed 
applicable to the known B[e] supergiants. 

The rotation velocities of B[e] supergiants. Our model is 
based on the assumption that B[e] supergiants are rapidly rotat- 
ing stars. What is the evidence for their rapid rotation? In fact 
not much is known about their rotation velocities. Due to their 
high density circumstellar medium, most of them do not show 
any photospheric absorption lines which might be used to de- 



rive a possible rotation speed. There are, however, three (out of 
15) B[e] supergiants in the Magellanic Clouds for which pho- 
tospheric absorption lines have been detected. From these line 
profiles only the projected stellar rotation (i.e. v sin z) can be de- 
rived with high accuracy. With their (often) poorly known incli- 
nations i, only lower limits to the real rotation speeds could be 
derived. These were found to be w > 0.35 and a> > 0.45 for the 
two LMC B[e] supergiants HenS93 and R66 (Gummersbach 
et al. [T995l Zickgraf 120061 respectively), and to - 0.8 for the 
SMC B[e] supergiant R 50 f Zickgraf l20TjQl Especially the lat- 
ter seems to rotate at a substantial amount of its critical velocity 
providing the basis for our research, although we cannot gen- 
eralize that all B[e] supergiants are rapidly rotating. 

Stellar luminosities and effective temperatures. A sum- 
mary of the stellar parameters (T e ff, L*, R*) of the MC B[e] 
supergiants is given e.g. in Zickgraf (2006 1. From this list it is 
obvious that the stellar luminosities and effective temperatures 
of the B[e] supergiants cover the range 10 4 5= log L/L Q ^ 10 6 
and 10000 < T eS [K] < 27 000, and our chosen values of the 
effective temperature fall well into this range while our lumi- 
nosities are taken as maximum values. The literature values for 
the B[e] supergiant effective temperatures, however, should be 
taken with caution because they have mainly been derived from 
fitting Kurucz model atmospheres to the observed spectral en- 
ergy distributions (see Zickgraf 1998 and references therein). 
However, these model atmospheres have been calculated un- 
der the assumption of spherically symmetric, non-rotating (i.e. 
uniformly bright) stars. If B[e] supergiants are indeed rapidly 
rotating, then for a proper determination of the mean (i.e. ob- 
servable) effective temperature a comparison with composite 
spectra should be undertaken (see e.g. Lovekin et al. 2006 1. Of 
course, to do so the inclination and the rotational velocities 
must be known which is usually not the case for B[e] super- 
giants. 

Mass loss rates. The models presented in Fig.|4] are for 
stars with polar mass fluxes between 1.0 x 10~ 5 and 2.0 x 
10~ 5 g s" 1 cm 2 . Therefore, the range in mass loss rates covered 
by these calculations extends from M m ; n = 3.4 x lO~ 5 M yr _1 
(critically rotating star with lower polar mass flux) to M max = 
1.3 x 10 _4 M o yr _1 (non-rotating star with higher polar mass 
flux), and is in good agreement with the known mass loss rates 
for MC B[e] supergiants, which range from about lO~ 5 M yr _1 
to about lO^Moyr" 1 (see Zickgraf 20©. 

The equatorial surface density. According to Eq. (TT2l . 
the surface density of a non-rotating star behaves as rin(R,) ~ 
M/(R 2 t v(R,)) = F m /v(R t ). For our model calculations we have 
chosen mass loss rates in the range reliable for B[e] super- 
giants, but we used some maximum values for the velocity 
(v{R„) = Voo) and stellar radius 2 . Therefore, for non-rotating 
stars the surface densities in our calculations provide some 
lower limits. 

For (especially rapidly) rotating stars, the situation be- 
comes more complicated. Now, the surface density is found 

2 There exists one exception in the literature: The (more or less) 
pole-on LMC star R66 is assumed to have a radius of 125 R Q (see 
Table 1 in Zickgraf 2006 1, while all other stars fall well below our 
adopted radius of 82 R Q . 
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to drop from pole to equator. On the other hand, B[e] super- 
giants are supposed to have circumstellar disks. These disks 
are much denser than what is found for the polar wind regions, 
with density contrasts on the order of p eq /p po i e - 100 .. . 1000. 
Models proposed to explain the formation of these high den- 
sity disks are the bi-stability mechanism (Lamers & Pauldrach 
119911 Pelupessy [2000 1 and wind compression (Bjorkman & 
Cassinelli ll993> . While the bi-stability mechanism in a rotating 
star might account for an increase by a factor of ~ 10 in equato- 
rial density only, the wind compression, especially for rapidly 
rotating stars (or more precisely the flow of material towards 
the equatorial plane), can be inhibited due to the appearance 
of a non-radial force provided by the radiation (Owocki et al. 
11996 ). Recently, the existence of a slow solution in line-driven 
winds of rapidly rotating stars has been found (Cure 2004 Cure 
et al. I2005> . Inclusion of the bi-stability jump resulted in an 
equatorial density enhancement (at least in the close vicinity 
of the star) by a factor of 100 - 1000, just what is needed to 
explain the disks of B[e] supergiants. However, these solutions 
have been found adopting a spherically symmetric star and ne- 
glecting gravity darkening, and it still needs to be confirmed 
that these slow solutions will also exist when gravity darkening 
is taken properly into account. 

Our models do not account for any density enhancements 
either due to bi-stability or due to wind compression. Such an 
increase in equatorial surface density by a factor 100 - 1000 
(for the same imput values) would result in a recombination 
distance even closer to the stellar surface, or would mean that 
the polar mass flux of the model star can be reduced by the 
same factor and the material would still recombine in the vicin- 
ity of the star. Such a lower mass flux (and hence a lower mass 
loss rate) might be desirable if the winds of B[e] supergiants 
are clumped. A clumpy wind in contrast to the assumed smooth 
density distribution is found to overestimate the mass loss rates 
(derived e.g. from Ha) by a factor of 10 or more (see e.g. Hillier 
120051 Bouret et al. l2005t . 

Even if the winds of B[e] supergiants will turn out to have 
lower mass loss rates, the ionization structure calculations pre- 
sented in this paper show that recombination of the equatorial 
wind material of rapidly rotating stars can take place at or at 
least close to the stellar surface. 

6. Conclusion 

We investigated the influence of rigid rotation on the surface 
and wind parameters of hot luminous stars, with emphasis on 
the non-spherical winds of B[e] supergiants. Since B[e] su- 
pergiants are known to have equatorial disks that show evi- 
dence for hydrogen neutral material in the vicinity of the stellar 
surface, the calculations are restricted to the equatorial plane. 
Due to the symmetric stellar radiation field (with respect to the 
equatorial plane) the problem of finding the recombination dis- 
tance reduces even to the 1 -dimensional case. The radiation 
field is however treated 2-dimensionally, to properly account 
for the latitude dependences of the parameters like effective 
temperature, wind velocity and density. The ionization balance 
equations are solved in a pure hydrogen plus helium wind. All 
assumptions made during our calculations have in common that 



the number of available ionizing photons at any location within 
the equatorial plane is overestimated. This results in a shift of 
the recombination distance to larger values, which means that 
we have calculated upper limits for the recombination distance. 

The major result is that despite the drop in equatorial sur- 
face density with increasing rotation velocity (neglecting any 
possible equatorial density enhancement due to bi-stability 
and/or wind compression), hydrogen recombines at (or close 
to) the stellar surface for rotating models with a polar mass 
flux F miPO ie ~ 1-5 x 10~ 5 g s^'cirT 2 and rotation velocities 
u> ^ 0.8 (see Fig.|4}. These mass fluxes correspond to mass 
loss rates M ^ 5 x 10~ 5 M o yr _1 for our chosen model stars 
with supergiant stellar and wind parameters. Since the mass 
loss rates for B[e] supergiants are found to lie in the range 
M = 10~ 5 . . . lO _4 M yr _1 we can expect that at least some 
of these stars might have hydrogen neutral equatorial material 
close to their stellar surface, given that they are indeed (rapidly) 
rotating stars, as is found for at least three of them. 
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